IASSNS-HEP-98/60 
PUPT-1802 
|hep-th/9806214 



Tachyons and Black Hole Horizons in Gauge 

Theory 

Daniel Kabat 1 and Gilad Lifschytz 2 

1 School of Natural Sciences 
Institute for Advanced Study 
Princeton, New Jersey 08540, U.S.A. 

kabat@ias . edu 

2 Department of Physics 
Joseph Henry Laboratories 
Princeton University 
Princeton, New Jersey 08544, U.S.A. 

gilad@pub.epl . princeton . edu 
Abstract 

Any probe which crosses the horizon of a black hole should be ab- 
sorbed. In M(atrix) theory, for 0-brane probes of Schwarzschild black 
holes, we argue that the relevant absorption mechanism is a tachyon 
instability which sets in at the horizon. We give qualitative arguments, 
and some quantitative large- ./V calculations, in support of this claim. 
The tachyon instability provides an attractive mechanism for infalling 
matter to be captured and thermalized by a Schwarzschild black hole. 
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1 Introduction 



It is natural to ask what M(atrix) theory |]J has to say about the physics 
of Schwarzschild black holes. Several groups have investigated this subject 
0, |, |, §, g, 0, g, 0,0,|Tl|,|g,g. By modeling the black hole as a Boltzmann 



gas of 0-branes, several key properties such as the entropy - mass relationship 
and the Hawking temperature have been derived, up to numerical coefficients 
of order unity. 

But the defining property of a black hole in classical gravity is the exis- 
tence of an event horizon. In this paper we address the question of how a 
horizon manifests itself in M(atrix) theory. One's first thought - that the 
size of the horizon is given by the range of eigenvalues of the matrices Xj - 
cannot be correct, since this range is expected to diverge in the large iV limit 
for all states in M(atrix) theory. 

In a scattering experiment the horizon area of a black hole can be mea- 
sured from its absorption cross section. In the context of M(atrix) theory 
this reflects the fact that a probe 0-brane should be absorbed by the black 
hole once it crosses the horizon. 

For Schwarzschild black holes in M(atrix) theory our proposal is that 
the mechanism responsible for this absorption is a tachyon instability: the 
probe 0-brane sees the horizon as a critical radius, at which strings stretched 
between the probe and the black hole become tachyonicQ This provides a 
precise prescription for identifying the horizon in M(atrix) theory. 

At a qualitative level, one can see that a state in gauge theory that 
gives rise to such a tachyon region will have the features expected of a black 
hole. Suppose a probe enters the tachyon region. Then the tachyons will 
start to condense, so the probe will be captured with very high probability. 
Moreover, the tachyonic modes couple to all the other degrees of freedom 
which make up the initial state, and the energy released as the tachyon 
condenses is available to be redistributed among all these degrees of freedom. 
Thus within the gauge theory we see that infalling matter gets absorbed and 
quickly thermalized when it enters a region with a tachyon instability. This 
is all directly analogous to the behavior expected of matter falling into a 
black hole. 

Although our treatment will focus on 11-dimensional Schwarzschild black 
holes, we expect the discussion of tachyon instability to apply more generally 



1 A horizon is a global property of spacetime, which cannot be detected by a local probe. 
More precisely our proposal is to associate the tachyon radius with a closed trapped surface, 
which we shall loosely refer to as a horizon. 
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to the horizon of non-extremal D-brane black holes, as well as to non-extremal 
black holes in the context of anti de Sitter space [|14[ . 



Extremal black holes must behave somewhat differently, however. A su- 
persymmetric probe of an extremal black hole should still be absorbed when 
it reaches the horizon, but the mechanism cannot be a tachyon instability. 
In section 3 we will give a brief discussion of string pair creation as a possible 
mechanism in the extremal case. 

An outline of this paper is as follows. In section 2 we review some facts 
about Schwarzschild black holes in M(atrix) theory and absorption mecha- 
nisms in D-brane scattering. In section 3 we argue that a tachyon instability, 
if present, would have the features expected of a Schwarzschild horizon in 
M(atrix) theory. In section 4 we study the tachyon instability for several ex- 
plicit M(atrix) configurations. Although we cannot directly study black hole 
states, the examples we consider are sufficient to show that tachyons are in- 
deed present in M(atrix) theory, in regions which can plausibly be associated 
with Schwarzschild horizons. Section 5 contains some conclusions. 



2 Preliminaries 

We review some relevant results concerning Schwarzschild black holes in 
M(atrix) theory and absorption mechanisms in D-brane scattering. 

2.1 Schwarzschild black holes in M(atrix) theory 

The most general object in M(atrix) theory is a collection of 0-branes with 
some strings stretching between them. The only truly stable states are the 
gravitons, but a black hole can be thought of as a long lived quasi-stable 
state. A Schwarzschild black hole has negative heat capacity, unlike any 
state in SYM. But in terms of light cone variables a black hole actually has 
a positive heat capacity if the number of non-compact dimensions is six or 
larger ||. This makes it possible to be described using SYM. 

Another feature of gravity that is relevant for black holes is that one 
can't put an arbitrary amount of energy in a bounded region. The upper 
limit for a given region is the black hole mass with the horizon area chosen 
appropriately. In M(atrix) theory we can see a similar effect. Fix the size 
of the matrices at iV x N, and choose some volume with typical length L. 
We want to see if there is a limit on the SYM energy, with the restriction 
that the object is contained within the designated volume. This is easy to 
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see at the classical level of the SYM. For the moment let us define the size 
to be at most the volume occupied by the eigenvalues of the matrices (which 
are not necessarily simultaneously diagonalizable). Now take nine diagonal 
matrices with some distribution of eigenvalues —L < A < L. This state has 
zero energy. The most general configuration with these eigenvalues can be 
obtained by acting with different unitary matrices on each of the original 
nine matrices. As the magnitude of each entry of a unitary matrix is smaller 
than 1, there is a limit on the entries of the resulting nine matrices. Then it 
is clear that the energy ~ Tr[Xj,X,,-] 2 is bounded. 

Currently the best understood model for Schwarzschild black holes is 
at the so-called BFKS point, that is when N ~ S |2|. At this point the 
temperature of the SYM is so low that quantum fluctuations are of the same 
order as thermal fluctuations. To go beyond this point, to N ^> S, requires 
a fully quantum mechanical treatment of the system. Following 0, one can 
model black holes with entropy ~ iV as follows. The basic ingredient is a 
state that has some expectation value for the matrices, Bi =< Xj > (this is 
referred to as the background). The SYM energy of this state is identified 
with the light-cone energy of the black hole, and the volume of the state is 
the volume of the black hole (we will come back to what we mean by volume 
of the state). The background and its fluctuations (which are modeled as a 
gas of N particles) contribute to the entropy and energy of the state. In fact 
most papers deal only with the thermodynamics of these fluctuations. 

For the entropy to come out correctly it is crucial that the statistics of 
the fluctuations be that of distinguishable particles. This issue is somewhat 
subtle. If all the background matrices commute (so that they can be simulta- 
neously diagonalized) then there is a residual gauge symmetry that permutes 
the diagonal elements, which makes the statistics that of identical particles. 
In this case the entropy is very small. To see this take a gas of free identical 
particles in d dimensions. The assumed kinematics is that each particle has 
momentum p±_ ~ 1/Rs set by the uncertainty principle and mass m = 1/R, 
where R$ is the Schwarzschild radius and R is the radius of the null circle. 
The volume of the gas is V = (Rs) d — N at the BFKS point. The energy of 
the gas is then E ~ NR/R 2 S . For a gas of N distinguishable particles in d 
space dimensions the entropy is 



This entropy is indeed ~ N. If the particles are identical then the volume V 
would be replaced with V/N and there would be very little entropy. 

If the background matrices do not commute (which means that there is 
a coherent state of strings stretched between different 0-branes) then we no 




(1) 
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longer have the residual gauge symmetry discussed above and thus do not 
have identical particles. The factor 1/N will not appear and the entropy 
will be order N. But the parameter which breaks the statistics symmetry is 
continuous, and it is not entirely clear at what value the effective statistics 
actually changes. For instance, if the energy in the fluctuations is much 
larger than the energy of the background then it is not clear why the above 
analysis should apply. 

Another issue is how to measure the size of a state. Clearly it cannot 
simply be the volume occupied by the 0-branes (although this will provide an 
upper bound) since, according to the holographic principle, this is expected 
to grow with N even for a graviton [0]. The most physical measure of the size 
of a state is through scattering experiments. For a black hole, in particular, 
it is natural to use the absorption cross-section as a definition of the horizon 
size. So before applying this definition to M(atrix) black holes we briefly 
review absorption mechanisms in D-brane scattering. 



2.2 Absorption in D-brane scattering 

In D-brane scattering there are two mechanisms which contribute to the 
imaginary part of the phase shift and lead to absorption. 

The first mechanism is string creation. When two D-branes have non-zero 
relative velocity, there is an amplitude to pair create two oppositely-oriented 
open strings stretched between them. Since velocity is T-dual to an electric 
field on the D-brane worldvolume, and open strings carry electric charge, this 



is dual to the phenomenon of pair production in an electric field [15| 



In the scattering amplitude, string creation arises as follows. The one- 
loop eikonal phase shift for DO - DO scattering has the form (b is the impact 
parameter and v is the relative velocity) |L6| 

*(&,„) ~ f^/w/M (2) 

Jo s sin sv 

At sv = nn the integrand has poles which contribute an imaginary part to the 
phase shift. The outgoing wavefunction of the scattered 0-brane is multiplied 
by e lS . Significant string creation is possible only for b < y/va', since for large 
b the imaginary part is exponentially suppressed, Im 5 ~ e~ b ^ va . This effect 
can be clearly seen in the SYM truncation of the one-loop amplitude, where 
the norm of the outgoing wavefunction is [T^, [T7| 



tmh8 Us ■ (3) 
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Of course for string creation to be possible there must be enough kinetic 
energy in the first place. So we see that there are two requirements for 
significant string creation to be possible, 

b < 2vW and b < 2na' — v 2 . (4) 

The second condition is the more stringent one when the relative velocity is 
small .0 

The second mechanism for absorption we refer to as tachyon instability 
Uy|. While there is always a term e~ b s in the integrand for the phase shift, 
it could happen that for large s the function f(s,v) behaves like e +c s . In 
eikonal scattering this happens when the oscillator ground state energy of 
the stretched string is negative. Then for b < c the proper time integral for 
the phase shift diverges. But by analytically continuing one sees that what 
actually happens is the phase shift acquires an imaginary part. 

Once a tachyon instability develops, the tachyon field will tend to con- 
dense. In a sense this condensation can be viewed as string creation, since 
a coherent state of stretched strings is produced (the strings just happen 
to have negative energy). But this picture of tachyons as stretched strings 
can be somewhat misleading in M(atrix) theory. In the context of SYM 
calculations the tachyon instability arises when the (mass) 2 matrix for the 
off-diagonal elements has a negative eigenvalue. Of course the full poten- 
tial Tr[Xj,Xj] 2 has a quartic piece that stabilizes the tachyon field at some 
finite value. When the tachyons condense, they lower the energy of the orig- 
inal configuration. They do this by making the off-diagonal elements larger, 
which naively means that additional strings are being added to the configu- 
ration. But after the tachyons condense, one can rediagonalize the matrices. 
One finds that the final configuration can be reinterpreted as having fewer 
stretched strings present than the initial state. Related phenomena have 
been discussed in string theory |L9|, ^0 . 



In view of this ambiguity, we will use the term 'string creation' only in 
reference to the pair production of open strings with positive energy arising 
from poles in the integrand. A point which will be important is that, even if 
no tachyon is present, string creation is enhanced when c 2 is positive. This is 
because the requirements for string creation (|J) are relaxed to have the form 



b < V^a'v + c 2 and b < 2na f I —v 2 + c ] (5) 



2 These are semiclassical estimates. The small velocity regime, where the second con- 
dition dominates, really calls for a fully quantum mechanical treatment of the DO - DO 
system. Such a treatment shouldn't change the qualitative behavior that string creation 
vanishes as the relative velocity goes to zero. We thank C. Bachas for raising this issue. 
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Thus string creation becomes possible for a wider range of impact parameters 
when c 2 is positive. 



3 Horizons in gauge theory 

An essential feature of classical black holes is that they have an event horizon. 
That is, there is a region of spacetime from which one cannot escape back to 
future infinity. This means that classically once an object enters the horizon 
it will be absorbed with certainty. 

We want to understand how this happens in gauge theory. In principle 
given the wave function of a black hole state in the gauge theory we could 
compute its absorption cross-section. This absorption could arise from one 
or both of the mechanisms discussed above. We first discuss Schwarzschild 
black holes, then briefly discuss the extremal case. 

3.1 Schwarzschild black holes 

To mimic the black hole, we expect absorption to set in at a relatively sharp 
radius, which should not depend on the nature of the probe. In particular 
the absorption cross section should be approximately independent of the 
incoming probe momentum. As we will see, this is a non-trivial requirement. 

We first consider string creation as a possible mechanism, and ask whether 
string creation by gas of 0-branes (treated as distinguishable particles) can 
be responsible for absorption by a Schwarzschild black hole at the BFKS 
point. The criterion is that if we send in a probe 0-brane with transverse 
momentum p±_ at an impact parameter smaller than the Schwarzschild radius 
it should be absorbed.^ This depends on the mean free path of a 0-brane in 
the gas of 0-branes making up the black hole. The mean free path is given 
by ' 

d = — (6) 
no- 
where n is the number density of the gas and a is the absorption cross-section. 
To estimate a we use the string creation conditions (f|). If the velocity is large 
then the absorption radius b ~ y 7 ^ and hence a ~ p\, while at small velocity 
b ~ v 2 so a ~ p l ®. The key point is that the probability of string creation 

3 At the BFKS point a single 0-brane can be used to probe a Schwarzschild black hole 
without much momentum transfer. 
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vanishes as the relative velocity goes to zero. 

The condition for absorption is that the mean free path be smaller than 
the Schwarzschild radius, d < Rs- Using the BFKS density n ~ 1/lp this 
translates into a lower bound on the transverse momentum of the probe. For 
large black holes it is the a ~ p 1 ^ estimate which controls the bound: a probe 
must have 

Px > (Rs)- 1/W (7) 

in Planck units to be absorbed. This is rather unsatisfactory. It would 
indicate that a Schwarzschild black hole with radius 10 38 / p ~ 1.6km can only 
absorb particles with momenta larger than 4 x 10~ 3 m p ~ 10 16 GeV. 

This seems to rule out pure string creation as a mechanism for absorption 
by the black hole. But one can consider the following hybrid mechanism. 
Recall from ([5]) that, even if no tachyon develops, string creation is enhanced 
when the oscillator ground state energy of a stretched string is negative. 
Basically this makes the string lighter so it is easier to pair produce them. 
One could imagine that this enhances string creation to the point where it 
is responsible for absorption by the black hole. 

Such an effect is very plausibly present in non-extremal M(atrix) black 
holes, as our examples in section 4 will show. But given that negative oscil- 
lator ground state energies are present and seem to play a significant role, it 
is natural to imagine that a full-blown tachyon instability might be present. 
So let us look more closely at the tachyon instability as an absorption mech- 
anism. 

Once a probe 0-brane enters the tachyon region, the tachyon mode will 
start to condense and roll down its potential. In the process it acquires 
kinetic energy^ But the tachyon field is coupled non-linearly to the other 
degrees of freedom making up the black hole via the Tr[X i; Xj] 2 interaction, 
so this energy is available for distribution among all black hole degrees of 
freedom. This is the process of absorbing and thermalizing the new matter 
that fell into the black hole. 

One might worry that the probe 0-brane will leave the unstable region 
before the tachyon condenses. The time scale for the tachyon to start rolling 
down the potential is t ~ V\/ — m tach y on- Unless the probe just grazes the 
horizon, in order not to be captured it must have a velocity of at least v ~ 
Rs/to. At such large velocities our quasistatic treatment of the background 

4 These estimates are for eleven dimensional black holes. In D spacetime dimensions a 
behaves as p^f 3 ^ 2 and p 2 ^ D 3 \ respectively. 

5 Quite unlike the string creation scenario where much of the initially available energy 
is used up in creating the string. 
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isn't reliable, although we note that at very high velocities string creation 
could be an important effect. One also might worry that the energy which is 
released by tachyon condensation could eject one of the other 0-branes from 
the black hole. But this is very unlikely at large N, when there are many 
degrees of freedom present that can share the energy. 

Of course 0-branes do escape from the black hole as Hawking particles 
. It would be interesting to understand how the tachyon instability affects 
this process. Suppose a particular 0-brane is liberated, in that the strings 
connecting it to the rest of the black hole happen to fluctuate to zero (either 
thermally or quantum mechanically). If this happens inside the tachyon 
region, then it is still classically almost impossible for the 0-brane to escape 
the black hole. It is as if the 0-branes have to tunnel out of the tachyon 
region in order to escape. 

Note that the tachyon instability does not depend sensitively on the rela- 
tive velocity of the probe and target - it is an effective absorption mechanism 
even if the probe has zero velocity. The tachyon instability does depend cru- 
cially on the structure of the target, however, and it is plausible that it could 
distinguish graviton states (which do not have horizons) from Schwarzschild 
black holes (which do), even though holography dictates that both gravitons 
and black holes in a sense grow with N. 

We feel that these properties make a tachyon instability an attractive 
mechanism for absorption by a Schwarzschild black hole. It seems natural 
to associate the tachyon region with the size of the black hole horizon. To 
make this plausible, we must show that tachyon instabilities can arise in 
macroscopically large regions in M(atrix) theory. We will do this in section 
4 through consideration of some explicit examples. 



3.2 Extremal black holes 

When a D-brane is used to probe an extremal black hole, absorption should 
still set in at the horizon. But the physical mechanism for absorption must be 
different from the Schwarzschild case, since a (supersymmetric) probe can't 
develop a tachyon instability. 



An explicit example of such a system was studied in ref. |2T|. These 
authors used a Dl-brane to probe an extremal black hole with three charges, 
Dl, D5, and momentum. One observation they made is that the horizon 
should be identified with the origin of the Coulomb branch of the probe 
gauge theory^. They studied geodesic motion on this moduli space, and 

6 This follows from the fact that in supergravity, in a coordinate system with the horizon 
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argued that certain geodesies which fall in towards p = are a signal of 
eventual absorption by the black hole. 

Unlike Schwarzschild black holes, the mechanism responsible for absorp- 
tion at the horizon of an extremal black hole cannot be a tachyon instability. 
Instead the mechanism is more plausibly one of string pair creation. Con- 
sider geodesic motion of a probe towards the origin of moduli space. The 
probe must have some non-zero velocity v. Then sufficiently close to the 
origin string creation will become possible and the moduli space approxima- 
tion will break downF]. String creation gives rise to an absorptive part of the 
scattering amplitude, and can lead to capture of the probe by the black hole 
(in gauge theory terms a transition from the Coulomb to the Higgs branch). 
As discussed above, string pair creation has a slow rate of thermalization, 
since no energy is released in the process. But this is consistent with the 
black hole being extremal. 



4 Examples of tachyon regions 

In this section we calculate the tachyon instability region (TIR) in some 
explicit M(atrix) configurations. Ideally we would perform this calculation 
for states corresponding to gravitons and Schwarzschild black holes. But 
our lack of precise knowledge of these wavefunctions forces us to turn to 
some simpler examples, which we nonetheless feel point out some general 
features of the TIR. Our purpose is to make it plausible that tachyons are 
present in M(atrix) theory, in macroscopic regions that can be associated 
with Schwarzschild black holes. 

The general calculational framework is as follows. We choose a.n N x N 
target configuration Bi, and place a 0-brane probe at a position foj. This 
corresponds to a background in M(atrix) theory 



< X, >= 



Bi 
h 



The covariant gauge- fixed action of the M(atrix) quantum mechanics is [23[ 
J dtTr (AX, AX, + ~[X i; X,] 2 + i* a D t *a - ^[X,, * 



9ym 

at p = 0, they found that the effective action for the probe has a term f(p) [dp 2 + p 2 dil 2 ) . 
But a term of the same form arises in the effective action on the Coulomb branch of the 
gauge theory, with p = interpreted as the origin of the Coulomb branch. Note the 
contrast to the non-extremal case, where the horizon is located away from the origin of 
moduli space |2^| . 

7 The one-loop estimate is that this occurs at p = y/va'. 
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"•"^gauge fixing ^"ghost • 

Expanding the action around the background to quadratic order gives the 
following (mass) 2 matrices for the off-diagonal degrees of freedom connecting 
the probe to the target (see for example |[27| ). 

E k (B k -h) 2 -WiBj-bj) 1 

2d t (B t - h) Sa E fc (B k - b k f + 2[B h B 3 ] \ 

^fermi = " ® 1 16xl6 + 1$ (A - 6,:) ® f + ^ [A, A] ® 7*' 

^host = E(A"^) 2 



We want to identify the probe positions which give rise to at least one 
negative (mass) 2 eigenvalue. Clearly M 2 host is positive definite so we do not 
need to worry about the ghosts. As long as the background is static M f 2 rmi 
is also the square of a Hermitian matrix and hence is positive definite. So 
tachyons do not arise in the fermion sector and we can ignore them as well.0 

This leaves the Yang-Mills degrees of freedom. If the background is time 
independent then Aq does not mix with Xj, and we only need to consider 
the following 9iV x 9N (mass) 2 matrix for the Xj degrees of freedom. 

(M 2 ) . . = Sij £ {B k - b k f + 2[B t , B 3 ] (9) 

k 

Even if the background does depend on time, dropping A in this way sets a 
lower bound on the size of the region with tachyons: given annxii Hermitian 
matrix with smallest eigenvalue A m j n , deleting one row and one column of the 
matrix produces a new (n — 1) x (n — 1) matrix whose smallest eigenvalue 
is larger than A min . So for the rest of the paper we will only analyze the 
mass matrix @. 

More generally, in situations where the target must be treated quantum 
mechanically, we can promote (§) to an operator equation. That is, we can 
promote the classical background Bi to a w(iV)-valued quantum operator Xj, 
and write f\ 

(M 2 ) . . = 6 {j J2 (X fc - b k f + 2[X, X,] (10) 

k 

8 In the context of string theory this is because the Ramond sector always has zero 
oscillator ground state energy. 

9 From now on Xj will denote the N x N operator corresponding to the target, not to 
be confused with the Ml (N + 1) X (N + 1) M(atrix) field. 
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where [•, •] denotes the commutator in u(N). A tachyon instability is present 
if the expectation value < i\)\M 2 \i\) > has a negative eigenvalue. This can be 
regarded as a proposal for defining the quantum horizon radius operator for 
a Schwarzschild M(atrix) black hole. 



4.1 An N = 2 example 

We start by considering the case N = 2. This is mainly a warm-up to illus- 
trate some properties of the TIR. An interesting feature which will emerge 
is that the TIR does not have to be connected. 

We consider backgrounds of the following form, involving two 0-branes 
and some stretched strings. 

Xx(0 = ai (t)a 3 (11) 
X 2 (i) = ^(ty (12) 

The a 1 are Pauli matrices. In general the operators ai(t), a 2 (t) obey Heisen- 
berg equations of motion, but for illustrative purposes we will assume that 
there is some classical time-independent background 

B x = ai a 3 (13) 
B 2 = a 2 u l (14) 

that approximates the properties of the state we are interested in. The 
constants ai, a 2 which appear in the background can be thought of as suitable 
time averages of quantum expectation values. 

The mass matrix for the bosons has the form 

M 2 = 

where C and D are 2x2 matrices. 

C = (a? + a^ + ^ + ^)l + 2aiM 3 + 2a 2 & 2 (7 1 (15) 
D = 2ia 1 a 2 a 2 (16) 

We want to find the region in (&i,fo 2 ) space where the probe has a tachyon 
instability, i.e. where M 2 has a negative eigenvalue. The boundary of this 
region can be found by looking for zeroes of the determinant. Without loss 
of generality we assume a 2 > a 2 , and find that the tachyon region is 




< b\ < \j\1a\al + 4(a| - a\)b\ -a\-b\ + a\ (17) 
max (o, a\ — a\ — Vl~2aia 2 ) <b\<a\ — a\ + Vl2aia 2 (18) 
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The maximum of b\ is 



max —= 2, \iYla\d2 — a 2 + a 1 \ (19) 



2 ' 

af 



If ai = a 2 then there is a tachyonic instability in a disk of radius b 2 = yYla\. 
If a 2 — a\ — \fY2a\<i2 > then there are two distinct regions in which there 
is a tachyonic instability. If a\ 3> a\ then the TIR consists of two regions, 
each roughly a circle of radius ~ Oi, centered at (0, a 2 ) and (0, —a 2 ). 

To look more precisely at the imaginary part, let us compute the phase 
shift of a 0-brane scattered from this configuration. To simplify things we 
will take the impact parameter b to be in the x 3 direction (and not as before 
in directions x 1 x 2 ), and give the probe 0-brane velocity v in direction x A . In 
this case there is still string creation and a tachyonic instability so this sim- 
ple example serves to illustrate the general case. Diagonalizing the (mass) 2 
matrices @ leads to the one-loop phase shift 

5(b,v) = j d ± e -^ 2 +a\+4)s_}_ ( 4 + 2 cogh 4 aia2S + 2cos2su (20) 
J s sin sv 

— 8 cos sv cosh 2aia 2 s) 

and there is a tachyonic instability for b 2 < 4aici2 — a 2 — a 2 ,. Notice that 
there is no tachyon in the fermionic sector (the last term in equation(pO|)). 
From equation (|20|) the norm of the outgoing wave function of the scattered 
0-brane can be computed to give 

(cosh 22f - cosh ^) 4 , 

VWt = \ 2 2 v ~^—o — 21 

W ' 8cosh 6 ^(cosh^ + cosh^) V ; 

where we have put r 2 = b 2 + a 2 + a\ and c = 2aia 2 . Let us look at equation 
(|2lD. If 02 = then we are reduced to scattering one 0-brane from two 
other 0-branes with the usual answer of tanh 8 {jib 2 /v). In the small v limit 
if r 2 < 4aia2 the norm of the wave function is ~ 0. Furthermore, even for 
small v and r 2 outside tachyonic region, if r 2 — 4aia2 ~ v then the norm of 
the wave function is ~ 1/4 signaling the enhancement of string creation close 
to (but outside) a TIR. 



4.2 Spherical membranes 

The next set of examples we consider are spherical membranes [^, . 
A round spherical membrane of radius r can be described in M(atrix) theory 
by setting 

B = \ 2 = 1,2,3 , 22 ^ 

) otherwise 
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The sphere is centered at the origin and extends in the directions x 1 x 2 x 3 . 
The matrices Jj are generators of SU(2) in the iV-dimensional representation, 
normalized to satisfy 

2 N 2 - 1 
[Ji, Jj] = i€ij k J k ^2 J i = — - — 1 . 

i 

We take the radius of the membrane to be much larger than the Planck 
length, so that this background can be treated semiclassically. The M(atrix) 
equations of motion imply that the membrane starts to collapse. Eventually 
it will fall through its Schwarzschild radius and form a black hole. But black 
hole formation seems out of reach of the semiclassical approximation, so we 
will study the membrane only when it is at its initial (maximal) radius, with 
^ = 0. This is feasible at large N because the time to collapse increases 
linearly with N (the formula is given in f27fl ). 



The M(atrix) background ( p2[) is a good approximation to a well-behaved 
state which is present in the full M-theory. It corresponds to a fundamental 
membrane, which has a degenerate horizon. In accord with this, we will 
verify that the region with tachyons for these matrices vanishes in the large 
N limit. This does provide evidence for the hypothesis that the TIR should 
be associated with the size of the horizon. 

We place a probe at a position 6j, which we take for simplicity to be non- 
zero only in the first three coordinates. Inserting the membrane background 
( |22|) into the (mass) 2 matrix (§), one finds that the non-trivial eigenvalues 
come from the 3iV x 3N block 

(M 2 ) % . = (r 2 (l - 1/iV 2 ) + |6| 2 ) 5 l3 l NxN - ^-rb ■ J5 tJ + i^e ijk J k 



where the indices i,j, k = 1, 2, 3 



Finding the exact eigenvalues of this matrix is a difficult problem. But, 
at leading order in 1/N, one can proceed as if the different SU(2) generators 
commute. Treating the Jj's as if they were numbers, and diagonalizing, one 
finds that the lowest eigenvalue comes from an N x N block given by 

A Q^2 



-^■J- W ^Jl + Jt + Jl 

In turn diagonalizing this block, one finds that the smallest eigenvalue of M 2 
is 

A min = (H-r(l-l/A0) 2 -2r 2 /iV 2 

where we dropped terms of order 1/iV 3 . So the smallest eigenvalue is negative 
if the probe is located in the region 

r(l-l/N)-y/2r/N < \b\ < r (1 - 1/iV) + y/2r/N . 
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That is, tachyons are present in a region which is centered at roughly the 
radius of the sphere, |6| ~ r, and extends in either direction to form a shell of 
thickness ~ r/N. In the large N limit the region with tachyons disappears. 
This is in accord with our expectations, since the fundamental membrane of 
M-theory has a degenerate horizon. 

The closely related example of a 0-brane probing a toroidal membrane 
was studied in p9l BUI . 



4.3 Gaussian Random Matrices 

To motivate this final example, we observe that a precise understanding of the 
relationship between M(atrix) theory tachyons and Schwarzschild black hole 
horizons would require finding the black hole wavefunction iphoieO^i, ■ ■ ■ > Xg). 
Then we could compute the expectation value of the (mass) 2 matrix 



<M 2 >= J d 9 X^ ole (X)M 2 Vvic(X) 



and find its lowest eigenvalue as a function of the position of the probe. If 
our proposed identification is correct, the lowest eigenvalue should become 
negative when the probe is exactly at the Schwarzschild radius. 

In attempting this, we are hampered by our poor understanding of the 
black hole wavefunction. Although many of its properties have been dis- 
cussed in the literature @, ^j, U, || [| [7|, ||, 0, [IT], 0, it is not understood 



precisely enough to make a calculation of < M 2 > possible. Nonetheless, 
several general properties which the wavefunction of a Schwarzschild black 
hole should satisfy are clear. 

1. At the BFKS point the wavefunction should have support on matri- 
ces Xj whose eigenvalues fall within a finite range. The range of the 
eigenvalues is presumably of order the Schwarzschild radius. 

2. A Schwarzschild black hole is spherically symmetric, so the wavefunc- 
tion should be invariant under SO(9) rotations Xj — > RijXj. 

3. The gauge symmetry of M(atrix) theory requires that the wavefunction 
be invariant under unitary transformations Xj — > ?7'Xj?7. 

Our strategy is to compute the distribution of eigenvalues of M 2 , not with 
respect to the true probability distribution given by solving M(atrix) theory 

P(X 1 ,...,X 9 )=^ le(X)^hole(X), 
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but rather with respect to the Gaussian probability distribution 



.P(Xi, . . . , Xg) = const, exp 




x 2a 2 

We regard the width of the Gaussian a as a freely adjustable parameter. 

The Gaussian ensemble in fact satisfies the general properties expected 
of a black hole listed above. The fact that the ensemble is 50(9) x U(N) 
invariant is obvious. Also, for large N the distribution of eigenvalues of any 
one of the matrices Xj is given by a Wigner semicircle 

1 



P{X) = 2^ VANa2 ~ ^ • (23) 

The large- N limit is taken with the 't Hooft coupling a 2 s = Na 2 held fixed, 
so that the eigenvalues of the matrices Xj all lie within the finite interval 

- 2a cff < A < 2a eS . (24) 



Since the Gaussian ensemble has some key features of a Schwarzschild 
black hole, we expect that studying the tachyon region in the Gaussian en- 
semble will point out some of the qualitative effects which would be present in 
a more realistic M(atrix) theory calculation. Indeed, an optimist could even 
hope to use the Gaussian ensemble as a reasonable variational approximation 



to the true black hole wavefunction at the BFKS point [31 



The practical advantage of Gaussian ensembles is that they are relatively 
easy to work with. Although we have not managed to find the exact region 
with tachyons in the Gaussian ensemble, we are able to put strict bounds on 
it as follows. Write the (mass) 2 matrix as a sum of two pieces. 

M 2 = A + B 

Aij = Sy ^2 X 2 , + 2[Xj, Xj] 



First we study the eigenvalue distribution of the matrix A, which corre- 
sponds to placing the probe at the origin (setting 6 = 0). It is convenient to 
work with dimensionless eigenvalues A defined by A = \/u 2 s - We introduce 
the resolvent 



a. 



2 

eff 



^( A ) = 7^ Tr ( 



9N a 2 s X-A' 
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Figure 1: The eigenvalue distribution pyi(A). Equivalently, the density of 
stretched string states vs. (mass) 2 for a probe at the center of a Gaussian 
ensemble. 

The normalized distribution of eigenvalues p A (A) can be obtained from the 
resolvent via 

Pa (A) = ImT^fA + ie) 

IT 

as a consequence of the identity = P.V.^ — Z7r5(a;). In the appendix 
we show that the resolvent, computed at leading order in l/N, satisfies the 
following cubic equation. 

73\K A + 584K 2 A + 9(9 - ~X)U A + 9 = (25) 

Of the three solutions to this equation, only one has a negative imaginary 
part. Solving for 1Z a and extracting the imaginary part gives the eigenvalue 
distribution shown in Fig. 1. As can be seen from the graph, the lowest 
eigenvalue of A is 

XT « -8.0a c 2 ff . (26) 
The key fact is that the lowest eigenvalue is negative. 

The distribution of eigenvalues of the matrix B can be obtained starting 
from the Wigner semicircle It is given by 



where the dimensionless probe position b is defined by b = \b\/a e s. Note that 
the eigenvalues of B lie in the interval 

b 2 - 46 < ~X B < b 2 + 4b . (27) 
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Having found the spectra of A and B separately, the lowest eigenvalue of 
M 2 = A + B can be bounded by 

\ min i \min ^ \min ^ \ min , \max /nn\ 
*A + A B S A M 2 < \ A + \ B ■ (28) 

If b = 0, so the probe is located at the origin, then the matrix B vanishes, 
and the eigenvalue distribution of M 2 is given in Fig. 1: a probe 0-brane 
placed at the center of a Gaussian ensemble has many tachyonic modes. 
As |6| increases the eigenvalues of M 2 tend to increase, and eventually a 
critical radius is reached at which there are no more tachyons. Applying 
the inequality ([28Q to the results (p6|) , ( pTj) we can say that the smallest 
eigenvalue of M 2 lies in the range 

b 2 - 46 - 8.0 < A™ 2 n < b 2 + Ab- 8.0 . (29) 

It then follows that the critical tachyon radius, at which the smallest eigen- 
value of M 2 vanishes, lies somewhere in the interval 1.47 < 6 cr i t < 5.47. 
Restoring units, this reads 

1.47<7eff < |&|crit < 5.47a eff . 

Recall from that the largest eigenvalue of the matrices Xj is 2a e g. 
So we see that the tachyon instability sets in at roughly the same radius as 
the size of the object defined by its range of eigenvalues. 

5 Conclusions 

We have argued that the horizon of a Schwarzschild black hole (or more 
precisely, a trapped surface) is seen by a 0-brane probe in M(atrix) theory 
as a critical radius at which a tachyon develops. This gives an interesting 
picture of the process by which such a black hole is formed. One starts with 
a collection of 0-branes located far from each other. With the right initial 
momenta they will approach each other and some string creation will take 
place. If enough strings are created then a tachyon region will form. If the 
size of the tachyon region is ^> l p then a macroscopic black hole has been 
created. If another 0-brane approaches the black hole, it will be absorbed 
when it enters the tachyon region. In the absorption process the tachyon 
rolls down its potential and, via its couplings to the other 0-branes making 
up the black hole, thermalizes the infalling matter. 

Since we do not have good control over the wavefunction of a Schwarzschild 
black hole, we cannot conclusively show that it has a tachyon region. But by 
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analyzing some examples where a tachyon instability arises in M(atrix) the- 
ory, we showed that the TIR has reasonable properties to be an indicator for 
a Schwarzschild horizon. We gave examples which have macroscopic tachyon 
regions, even in the large N limit. We also studied membrane states, where 
the TIR vanishes at large N. 

The most interesting example we analyzed was the tachyon instability 
for Gaussian random matrices. In this large- N calculation we showed that 
the size of the tachyon region is roughly the same as the spread of 0-brane 
positions. It is conceivable that the Gaussian ensemble is a good variational 
approximation to the true wavefunction of a Schwarzschild black hole at the 
BFKS point. At least in the non-supersymmetric case, Gaussians have been 
shown to give a fairly good approximation for large- iV Yang-Mills quantum 
mechanics [3~2"||. The crucial test of this idea is to include the fermions and 



compute the mass of a typical state in the Gaussian ensemble [31 



We expect that much of our discussion of tachyon instabilities can be car- 
ried over to other non-extremal black holes built out of D-branes. And given 
the conjectured relation between finite temperature SYM and black holes 
in AdS space [Q, our proposal suggests that the phenomenon of tachyon 
instability can also be used as a gauge theory indicator for the horizon of a 
non-extremal AdS black hole. 
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A Gaussian Random Matrices 

In this appendix we review some techniques of random matrix theory, and 
derive the cubic equation (|25D satisfied by the resolvent 71 a- 



The standard reference on random matrix theory is the book by Mehta 
3^|. The methods we will use were explained to us by Misha Fogler. Al- 



though we are not sure of original references similar methods appear in the 
literature in p4|. 



To introduce the method, we start by considering a single N x N Hermi- 
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tian matrix X chosen at random from the probability distribution 

P(X) = const, exp (-^Tr(X 2 )) . 
We are interested in finding the ensemble average of the Green's function 



G(A) ' 



A-X 



G(A) can be expanded in perturbation theory, in terms of a 'propagator' 
Go = t1 and an 'interaction' X, as 

G = Go ~\~ GqX-Gq -\- GoXGoXGo ~\~ ■ ■ ■ • 

To compute < G > we must contract the X matrices in all possible ways 
using the two-point function 

< X a ;,X cc i >= a 2 5 a dSbc ■ 



We take the large N limit holding the 't Hooft coupling a 2 s = Na 2 fixed. 
Then the diagrams which contribute to < G > at leading order in 1/N are 
the rainbow diagrams, in which X propagators do not cross. The same set 



of diagrams appear in the solution of large- iV QCD in 1+1 dimensions |p5 
They can be summed by setting 

< G >= G + G S < G > 

where the amputated, 1PI (with respect to G ) self-energy satisfies 



X < G > X 

LiVxiv 



a 2 t NxN Ti <G> 



Combining these two equations yields a quadratic equation satisfied by the 
resolvent 1Z = ^Tr < G >. Working in terms of the dimensionless quantity 
A = A/cr e ff the resolvent satisfies 

U 2 - ~\K + 1 = . 

Taking the appropriate solution to this equation, one finds that 7£(A) has 
a negative imaginary part for — 2 < A < 2, and that the distribution of 
eigenvalues of X is 

1_ _,r . 1 



p(A) = — Imft(A + ie) = — V4 - A 2 . 

7T 2,71 

This is Wigner's semicircle distribution (p3|), written in terms of A. 
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Next we consider a slightly more involved example, which has many of 
the features of the full problem. Let X and Y be two independent Hermi- 
tian matrices chosen at random from identical Gaussian distributions. To 
determine the eigenvalues of the commutator i[X, Y] we introduce as before 

G < A > = x^kr] 

< G > = G + G £ < G > 



The 1PI self-energy is given by 



E = X < YGY > X + Y < XGX > Y 

= a 2 t NxN Ti < XGX > +a 2 l 7Vx 7vTr < YGY > 

The quantities appearing on the right hand side in turn satisfy 



< XGX > = X<G>X + X<G>X< YGY > X < G > X 
= <r 2 l NxN Tr < G > +cr 4 < YGY > (Tr < G >) 2 



< YGY > = Y < G > Y + Y < G > Y < XGX > Y < G > Y 

= a 2 l NxN Tr <G>W< XGX > (Tr < G >f 

Other terms which could appear on the right hand side either vanish iden- 
tically or are sub leading in 1/N. This system of equations can be combined 
into a single cubic equation for TZ = ^Tr < G >. In terms of A = A/cr 2 ff , it 
reads 

Aft 3 + TZ 2 - ~\TZ + 1 = . 

Choosing the appropriate solution to this equation, one finds that the eigen- 
value distribution of z[X, Y] is symmetric about zero, and is non- vanishing 



only for A between ±| \J22 + 10\/5. 

Finally, we turn to the spectrum of the matrix 

A a = 5 ij Yl x fc + 2 t Xi > x j] 

k 

where the Xj are independent, identically distributed N x N Gaussian ran- 
dom matrices. We let i,j — 1, . . . , d so that A is an Nd x Nd matrix. As 
above we set 

1 



G(A) 



X-A 

< G > = G + G S < G > 



10 Note that compared to the previous example an extra power of a e g is present in the 
definitions of TZ and A. 
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where now Go = t1 NdxNd- The self-energy is given by 



where the contractions in the second term mean that the two outermost Xj's 
which appear in A ik , Aij are to be contracted. This equation can be simplified 
somewhat by making use of SO(d) invariance, which implies that we can set 

G{j ~5ijG ' —dij^jfck . 

Then one has 



— X fc X k 5ij + — < A ik GuA k j > 



gd — 8 

Nda 2 5ijl NxN H — 2 — cr 2 5 i: ,lAr X ArTr < X fc G w X fc > 



The quantity appearing on the right hand side in turn satisfies 



< X/jGz/Xjfc > _ Xfc < Gu > Xfc + Xjt < Gi m > A 

mil < G n i > X 



k 



+ Xfc < Gi m >< A mn G no A op >< G p i > X/j 

dcr 2 ljvxivTr < G k k > +ct*1nxn (Tr < G k k >) 
9d-8 



-a < XfeGiiXfc > (Tr < G mm >] 



There is a unique choice of contractions which contribute at leading order in 
1/N, since the X propagators cannot cross and the expectation value of an 
odd number of Xj's vanishes. The second equality is most easily obtained 
with the help of SO(d) invariance. 

This set of equations can be combined into a single cubic equation satisfied 

2 ~ 

by the resolvent TZa = ^ff Tr < G >. In terms of A = the equation reads 

efT 

(9d - 8)\K 3 A + (9d - 8)(d - l)U 2 A + d(d - \)TZ A + d = 0. 
Setting d = 9, this is the result (E5|) given in the text. 
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